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/-y^ ' A formalism to evaluate the Sivers function, developed for calculations in con- 

^~^ ' stituent quark models, is applied to the Isgur-Karl model. A non-vanishing 

Sivers asymmetry, with opposite signs for the u and d flavor, is found; the 
Burkardt sum rule is fulfilled up to 2%. Nuclear effects in the extraction of 
I J neutron single spin asymmetries in semi-inclusive deep inelastic scattering off 

|^N( ' ^He are also evaluated. In the kinematics of JLab, it is found that the nuclear 

jH , effects described by an Impulse Approximation approach are under control. 

Keywords: DIS, transversity, neutron structure. 



1. The Sivers function in Constituent Quark Models 

The partonic structure of transversely polarized nucleons is still an open 
problem.^ Semi-inclusive deep inelastic scattering (SIDIS) is one of the pro- 
posed processes to access the parton distributions (PDs) of transversely po- 
larized hadrons. SIDIS of unpolarized electrons off a transversely polarized 
target shows "single spin asymmetries" (SSAs),'^ due to two physical mech- 
anisms, whose contributions can be distinguished,'^"^ i.e. the Collins^ and 
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Fig. 1. The contributions to the Sivers function in the present approach. 



the Sivers^ mechanisms. The form.er is due to parton final state interactions 
(FSI) in the production of a hadron by a transversely polarized quark. The 
Sivers mechanism leads to a SSA which is the product of the unpolarized 
fragmentation function with the Sivers PD. The latter describes the num- 
ber density of unpolarized quarks in a transversely polarized target: it is a 
time-reversal odd, Transverse Momentum Dependent (TMD) PD. From the 
existence of leading-twist Final State Interactions (FSI) J'^ a non- vanishing 
Sivers function has been explained as generated by the gauge link in the 
definition of TMDs,^'^" whose contribution does not vanish in the light- 
cone gauge, as happens for the standard PD functions. Recently, the first 
data of SIDIS off transversely polarized targets have been published, for 
the proton^^ and the deuteron.^^ It has been found that, while the Sivers 
effect is sizable for the proton, it becomes negligible for the deuteron, so 
that apparently the neutron contribution cancels the proton one, showing 
a strong flavor dependence of the mechanism. Different parameterizations 
of the available SIDIS data have been published, ^^"^^ still with large error 
bars. Since a calculation from first principles in QCD is not yet possible, 
several model evaluations have been performed, e.g. in a quark-diquark 
model;^'^'^^ in the MIT bag model;^^ in a light-cone model;^^ in a nuclear 
framework, relevant to proton-proton collisions. ^^ We here describe a Con- 
stituent Quark Model (CQM) calculation of the Sivers function.^" CQM 
calculations of PDs are based on a two steps procedure. ^^ First, the ma- 
trix element of the proper operator is evaluated using the wave functions 
of the model; then, a low momentum scale, /Iq, is ascribed to the model 
calculation and QCD evolution is used to evolve the observable calculated 
in this low energy scale to the scale of DIS experiments. Such procedure 
has proven successful in describing the gross features of PDs^^and GPDs,^'^ 
by using different CQMs, e.g. the Isgur-Karl (IK) model. ^* Besides the fact 
that it successfully reproduces the low-energy properties of the nucleon. 
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the IK model contains the one-gluon-exchange (OGE) mechanism.^"'' In the 
present calculation, with respect to calculations of PDs and GPDs, the lead- 
ing twist contribution to the FSI has to be taken into account. The main 
approximations have been: i) only the valence quark sector is investigated; 
ii) the leading twist FSI are taken into account at leading, OGE, order, 
which is natural in the IK model; iii) the resulting interaction has been 
obtained through a non-relativistic (NR) reduction of the relevant opera- 
tor, according to the philosophy of constituent quark models,^''' leading to a 
potential Vnr. The Sivers function for a proton polarized along the y axis 
and for the quark of flavor Q, f:^{x, kx), takes the form (cf. Fig. 1 for the 
labels of the momenta and hclicities): 

= 3 l-ig^^ J dhdh^5{k+ - xP+)S{hT + qr - M^I^U) 
where g is the strong coupling constant, M the proton mass, and 
X«W ^ J2 ^lfs^^^(k3,m3;ki,mi;P-k3-ki,n 

nil ,m^,m3,m3 

X VNR{ki,k3,ci) 

X ^sf,s,=~i[h + q,m'^; ki- q,m[; P ~k2~ki,mn) . (2) 



Using the spin-flavor wave function of the proton in momentum space, 
$s/, corresponding to a given CQM, the Sivers function, Eq. (1), can be 
evaluated. From Eq. (2), one notices that the helicity conserving part of 
the global interaction does not contribute to the Sivers function. Besides, 
in an extreme NR limit, it turns out to be identically zero: in our scheme, 
it is precisely the interference of the small and large components in the 
four-spinors of the free quark states which leads to a non-vanishing Sivers 
function. This holds even from the component with ^ = of the target wave 
function. While, in other approaches, ^^ these interference terms arise due 
to the wave function, they are produced here by the interaction. 

The above-described formalism is now applied to the IK model. The 
detailed procedure and the final expressions of the Sivers function in this 
model can be found in Ref.^° To evaluate numerically Eq. (1), g (i.e. as(Q^)) 
has to be fixed. The prescription^^ is used to fix ^q, according to the amount 
of momentum carried by the valence quarks in the model. Here, assuming 
that all the gluons and sea pairs in the proton are produced perturbatively 
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Fig. 2. 



Left (right): the quantity f-^j; '^^ (^)i Eq. (3). Dashed eurve: IK at /^jj. FuU 



eurve: the evolved distribution at NLO. Patterned area: parameterization by (see text). 



according to NLO evolution equations, in order to have ~ 55% of the mo- 
mentum carried by the valence quarks at a scale of 0.34 GeV^ one finds that 
^il ~ 0.1 GeV2 if A^^g ~ 0.24 GeV. This yields a,(^g)/(47r) ~ O.lS.^i The 
results of the present approach for the first moments of the Sivers function, 
defined as 



/: 



-(1)2^^^ - 



IT 



(x) = d^h 



'2M2 



f^^{x,kT) , 



(3) 



are given by the dashed curves in Fig. 2. They are compared with a param- 
eterization of the HERMES data, taken at Q^ = 2.5 GeV^ : The patterned 
area represents the 1 — cr range of the best fit proposed in Ref.^'* The magni- 
tude of the results is close to that of the data, although they have a different 
shape: the maximum (minimum) is predicted at larger values of x. Actu- 
ally ^Q is much lower, Q^ = 2.5 GeV^. A proper comparison requires QCD 
evolution of TMDPDs, what is, to large extent, unknown. We nevertheless 
perform a NLO evolution of the model results assuming, for /-^gi {x), the 
same anomalous dimensions of the unpolarized PDFs. From the final result 
(full curve in Fig. 2), one can see that the agreement with data improves 
dramatically and the trend is reasonably reproduced at least for x > 0.2. 
Although the performed evolution is not exact, the procedure highlights 
the necessity of evolving the model results to the experiment scale and it 
suggests that the present results could be consistent with data, still affected 
by large errors. 

Properties of the Sivers function can be inferred from general principles. 
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The Burkardt Sum Rule (BSR)^^ states that, for a proton polarized in the 
positive y direction, X]q=m di^'x) — ^i^h 



j\x j dkr^ft^^ixM) . (4) 



Jo 

and must be satisfied at any scale. Within our scheme, at the scale of the 
model, it is found {k^) == fO.85 MeV, {kf) = -11.25 MeV and, in order 
to have an estimate of the quality of the agreement of our results with 
the sum rule, we define the ratio r — |(fcf) + {k^)\/\{k'^) — (fc")| obtaining 
r ~ 0.02, so that we can say that our calculation fulfills the BSR to a 
precision of a few percent. One should notice that the agreement which is 
found is better than that found in other model calculations, ^^'^^ especially 
for what concerns the fulfillment of the Burkardt Sum Rule. 

2. The Sivers function from neutron (^He) targets 

As explained in the previous section, the experimental scenario which arises 
from the analysis of SIDIS off transversely polarized proton and deuteron 
targets^^'^^ is puzzling. The data show an unexpected fiavor dependence 
in the azimuthal distribution of the produced pious. With the aim at ex- 
tracting the neutron information to shed some light on the problem, a 
measurement of SIDIS off transversely polarized ■^He has been addressed, ^^ 
and two experiments, aimed at measuring the azimuthal asymmetries in the 
production of leading tt* from transversely polarized '^He, are forth-coming 
at JLab.^® Here, a realistic analysis of SIDIS off transversely polarized "^He^^ 
is described. The expressions of the Collins and Sivers contributions to the 
azimuthal Single Spin Asymmetry (SSA) for the production of leading pi- 
ous have been derived, in impulse approximation (lA), including the initial 
transverse momentum of the struck quark. The final equations are involved 
and they are not reported here. They can be found in.^^ The same quantities 
have been then evaluated in the kinematics of the JLab experiments. Wave 
functions'^° obtained within the AV18 interaction'^^ have been used for a re- 
alistic description of the nuclear dynamics, using overlap integrals evaluated 
in Ref.,^^ and the nucleon structure has been described by parameteriza- 
tions of data or model calculations. ^•^■'^^ The crucial issue of extracting the 
neutron information from ^He data will be now discussed. As a matter of 
facts, a model independent procedure, based on the realistic evaluation of 
the proton and neutron effective polarizations in "^He,'^^ called respectively 
Pp and Pn in the following, is widely used in DIS to take into account ef- 
fectively the momentum and energy distributions of the bound nucleons in 
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^He. It is found that the same extraction technique can be apphed also in 
the kinematics of the proposed experiments, although fragmentation func- 
tions, not only parton distributions, are involved, as it can be seen in Figs. 
1 and 2. In these figures, the free neutron asymmetry used as a model in 
the calculation, given by a full line, is compared with two other quantities. 
One is: 

1 
d„ 



Al, ^ ^Ar^'^ , (5) 



A„ 



where i stands for "Collins" or "Sivers", A'^^^''^ is the result of the full 
calculation, simulating data, and dn is the neutron dilution factor. The 
latter quantity is defined as follows, for a neutron n (proton p) in ^He: 

E, e^/«'"(P) (x) £'«''* (z) 

and, depending on the standard parton distributions, /' (x), and frag- 
mentation functions, D'^'^ (z), is experimentally known (see^^ for details). 
j4^ is given by the dotted curve in the figures. The third curve, the dashed 
one, is given by 

-l-(Ar--2p,d,A;-^-) , (7) 

i.e. "^He is treated as a nucleus where the effects of its spin structure, of 
Fermi motion and binding, can be taken care of by parameterizing pp and 
Pn- One should realize that Eq. (5) is the relation which should hold between 
the "^He and the neutron SSAs if there were no nuclear effects, i.e. the "^He 
nucleus were a system of free nucleons in a pure S wave. In fact, Eq. (5) 
can be obtained from Eq. (7) by imposing pn = 1 and pp = 0. It is clear 
from the figures that the difference between the full and dotted curves, 
showing the amount of nuclear effects, is sizable, being around 10 - 15 % 
for any experimentally relevant x and z, while the difference between the 
dashed and full curves reduces drastically to a few percent, showing that 
the extraction scheme Eq. (7) takes safely into account the spin structure of 
^He, Fermi motion and binding effects. This important result is due to the 
kinematics of the JLab experiments, which helps in two ways. First of all, 
to favor pions from current fragmentation, z has been chosen in the range 
0.45 < z < 0.6, which means that only high-energy pions are observed. 
Secondly, the pions are detected in a narrow cone around the direction of 
the momentum transfer. As it is explained in,^^ this makes nuclear effects 
in the fragmentation functions rather small. The leading nuclear effects 
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Fig. 3. Left (right) The model neutron Collins (Sivcrs) asymmetry for tt^ production 
(full) in JLab kinematics, and the one extracted from the full calculation taking into 
account the pp (dashed), or neglecting it (dotted). The results are shown for z=OAb and 
Q^ = 2.2 GcV^, typical values in the kinematics of the JLab experiments. 

are then the ones affecting the parton distributions, already found in DIS, 
and can be taken into account in the usual way, i.e., using Eq. (7) for the 
extraction of the neutron information. In the figures, one should not take the 
shape and size of the asymmetries seriously, being the obtained quantities 
strongly dependent on the models chosen for the unknown distributions.^^ 
One should instead consider the difference between the curves, a model 
independent feature which is the most relevant outcome of the present 
investigation. Eq. (7) is therefore a valuable tool for the experiments.^^ 
The evaluation of final state interactions effects and the inclusion of more 
realistic models of the nucleon structure are in progress. 
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